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In this paper, we study the discrete-time approximation of mul- 
tidimensional reflected BSDEs of the type of those presented by Hu 
and Tang [Probab. Theory Related Fields 147 (2010) 89-121] and 
generalized by Hamadene and Zhang [Stochastic Process. Appl. 120 
(2010) 403-426]. In comparison to the penalizing approach followed 
by Hamadene and Jeanblanc [Math. Oper. Res. 32 (2007) 182-192] 
or Elie and Kharroubi [Statist. Probab. Lett. 80 (2010) 1388-1396], 
we study a more natural scheme based on oblique projections. We 
provide a control on the error of the algorithm by introducing and 
studying the notion of multidimensional discretely reflected BSDE. 
In the particular case where the driver does not depend on the vari- 
able Z, the error on the grid points is of order i — e, e > 0. 



1. Introduction. The main motivation of this paper is the discrete-time 
approximation of the following system of reflected backward stochastic dif- 
ferential equations (BSDEs) 

'Y; = g^XT)+ r f{x,x,zi)ds- r ZldW. 



(1.1) 



y/>max{y/-c^^(X,)} 

y/-max{y/ -c^^(Xi)} 



dKl = 0, 



o<t<r, 
o<t<r, 



Received April 2010; revised February 2011. 

^Supported by the "Chaire Risque de credit," Federation Bancaire Frangaise. 
^Supported by the the Finance and Sustainable Development Chair sponsored by EDF, 
CACIB and CDC. 

AMS 2000 subject classifications. Primary 93E20, 65C99; secondary 60H30. 
Key words and phrases. BSDE with oblique reflections, discrete time approximation, 
switching problems. 

This is an electronic reprint of the original article published by the 
Institute of Mathematical Statistics in The Annals of Applied Probability, 
2012, Vol. 22, No. 3, 971-1007. This reprint differs from the original in pagination 
and typographic detail. 



1 



2 



J.-F. CHASSAGNEUX, R. ELIE AND I. KHARROUBI 



where X := /, g and are Lipschitz functions and X is 

the solution of a forward stochastic differential equation (SDE). 

These equations are linked to the solutions of optimal switching problems, 
arising, for example, in real option pricing. In the particular case where / 
does not depend on {Y,Z), a first study of these equations was made by 
Hamadene and Jeanblanc [13]. They derive existence and uniqueness of so- 
lution to this problem in dimension 2. The extension of this result to optimal 
switching problems in higher dimension is studied by Djehiche, Hamadene 
and Popier [9], Carmona and Ludkovski [5], Porchet, Touzi and Warin [23] 
or Pham, Ly Vath and Zhou [22] for an infinite time horizon consideration. 
In this last paper, the resolution of optimal switching problems relies mostly 
on their link with systems of variational inequalities. 

Considering deterministic costs, Hu and Tang [15] derive existence and 
uniqueness of solution to this type of BSDE and relate it to optimal switch- 
ing problems between one-dimensional BSDEs. Extensions developed in [14] 
and [8] cover, in particular, the existence of a unique solution to the BSDE (1.1). 
Recently two of the authors related in [11] the solution of (1.1) to correspond- 
ing constrained BSDEs with jumps. As presented in [12], this type of BSDE 
can be numerically approximated combining a penalization procedure with 
the use of the backward scheme for BSDEs with jumps; see [2]. Unfortu- 
nately, no convergence rate is available for this algorithm. We present here 
a more natural discretization scheme based on a geometric approach. For 
any t <T, all the components of the Yt process are interconnected, so that 
the vector Yj lies in a random closed convex set Q{Xt) characterized by the 
cost functions {d^)ij^x- The vector process Y is thus obliquely reflected on 
the boundaries of the domain Q{X) and we approximate these continuous 
reflections numerically. 

As in [1, 6, 18], we first introduce a discretely refiected version of (1.1), 
where the refiection occurs only on a deterministic grid K = {rg := 0, . . . , 
r« := T}:Y^ = Y^ := giXr) G Q{Xt), and, for j < k - 1 and t e [r^, r^+i), 

\ Y^ = + p(Xi, y,«)i|,es}, 

where V{Xt,-) is the oblique projection operator on Q{Xt), for t<T. Ex- 
tending the approach of Hu and Tang [15], we observe that the solution 
to (1.2) interprets as the value process of a one-dimensional optimal BSDE 
switching problem with switching times belonging to K. This allows us to 
prove a key stability result for this equation. We control the distance be- 
tween (y^, Z^) and (Y, Z) in terms of the mesh of the reflection grid. Due 
to the obliqueness of the reflections, the direct argumentation of [1, 6] does 
not apply. Using the reinterpretation in terms of switching BSDEs, we flrst 
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prove that approaches Y on the grid points with a convergence rate 
of order ^ — £, e > uniformly in 9?, whenever the cost function is Lips- 
chitz and / is bounded in z (see Theorem 5.2). Imposing more regularity 
on the cost functions, we control the convergence rate of {Y^, Zf^)(j<t<T to 
{Yt,Zt)o<t<T (see Theorem 5.3). 

We then consider a Euler type approximation scheme associated to the 
BSDE (1.2) defined on vr = {to, ...,*„} by := £f(XJ) and, for i € {n - 
1,...,0}, 

- u)-'E[Y^;:iWu_^, - WuY \ ^ui 

i+1 



(1.3) 



z 



Y, 



»,71 



J•iJ + (^,+l-^,)/(x,-,y,f^zJ' 



where X'^ is the Euler scheme associated to X. It is now well known (see, 
e.g., [4, 24]), that the convergence rate of the scheme (1.3) to the solution 
of (1.2) is controlled by the regularity of {Y,Z) through the quantities 



E 



E 



\Y, 



Y,^\^dt 



and E 



i<n 



\zl 



ZffAt 



-E[//^+^Zf dt| JiJ, for i< 



n. 



with Zf^ = T- .^11. 

Using classical Malliavin differentiation tools, we prove a representation 
for Z^, extending the results of [1, 6] to the system of discretely reflected 
BSDEs (1.2). We deduce the expected regularity results on {Y^,Z^) and, 
using the techniques of [7], Chapter 3, we obtain in a very general setting 
the convergence of (1.3) to (1.2). However, due to the obliqueness of the 
reflections, the projection operator V{X, •) is only L-p-Lipschitz with L-p := 
\Q > 1, leading to a convergence rate controlled by |Lp|'^(|7r|^/^ + K^/^|7r|"^/^), 
where we recall that k is the number of points in the reflection grid K. The 
term |ivp|'^ can be very large even for small k. and leads to a poor logarithmic 
convergence rate when passing to the limit k — t- oo for the approximation 
of (1.1). In the particular case where / does not depend on z, we are able 
to get rid of the \L'p\'^ term. 

Our innovative approach relies on the use of comparison results to get 
a control of the involved quantities: 

• we interpret the solution of (1.2) as a value process of an optimization 
problem, which allows us to get a control of the distance between the 
continuously and discretely reflected BSDEs; 

• we introduce a convenient auxiliary process dominating both solutions (1.2) 
and (1.3) to get a control of the distance between these quantities. 

Combining the previous estimates, we deduce the convergence of the discrete 
time scheme (1.3) to the solution of (1.1) with a convergence rate of order 
^ — £, e > 0, on the grid points, whenever K = vr and / is independent of Z. 
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Whenever the cost functions are constant, ah the previous estimates hold 
true with e = 0. We want to emphasize that ah these resuhs are obtained 
without any assumption on the nondegeneracy of the volatihty matrix a. 

The rest of the paper is organized as fohows. In Section 2, we introduce 
the notion of discretely obliquely reflected BSDEs, connect it with opti- 
mal switching problems and give the fundamental stability result. Section 3 
focuses on the regularity of the solution to this new type of BSDE. This 
analysis leads to precious estimates allowing us to deduce the convergence 
of the associated discrete time scheme (see Section 4). Afterward, Section 5 
focuses on the extension to the continuously reflected case and provides 
a convergence rate of the discretely reflected BSDE to the continuously one, 
whenever the driver / is bounded in the variable Z. The global error of the 
scheme is provided at the end of this section. Some a priori estimates are 
reported in the Appendix. 

Notation. Throughout this paper we are given a finite time horizon T 
and a probability space {VL,F,W) endowed with a d-dimensional standard 
Brownian motion W = {Wt)t>o- The filtration F = [Ft)t<T generated by 
the Brownian motion is supposed to satisfy the usual conditions. Here, *P 
denotes the cr-algebra on [0, T] x generated by F-progressively measurable 
processes. Any element a; G with ^ G N will be identified to a column 
vector with ith component x* and Euclidean norm |a::|. For x,y G M^, x • y 
denotes the scalar product of x and y, and x' denotes the transpose of x. 
We denote by ^ the component by component partial ordering relation on 
vectors. A^"^'"^ denotes the set of real matrices with m lines and d columns. 
We denote by the set of functions from to M with continuous and 
bounded derivatives up to order k. For a function f £ C^, Vxf denotes 
the Jacobian matrix of / with respect to x. For ease of notation, we will 
sometimes write Ef[-] instead of E[-|J^j], t G [0,r]. In the following, we shall 
use the notation without specifying the dimension nor the dependence in 
u £ Q when it is clearly given by the context. Finally, for any p > 1, we 
introduce the following: 

• the set of real-valued cddldg^ 'P-measurable processes Y = (yi)o<t<T 
satisfying : = E[supo<t<T < oo. 

• the set T-L^ of M'^-valued ^-measurable processes Z = {Zt)o<t<T such that 
\\Z\\^,:=E[{J^\Zt\'dtr/Y/P<oo. 

• the closed subset of consisting of nondecreasing processes K satis- 
fying Ko = 0. 

In the sequel we denote by a constant whose value may change 
from line to line but which depends only on L. We use the notation C£ 
whenever it depends on some other parameter p > 0. 



^French acronym meaning right continuous with left limit. 
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2. Discretely obliquely reflected BSDE. In the beginning of this section 
we define and study discretely obliquely reflected BSDEs in a general set- 
ting. In particular, we show how their solutions relate to the solutions of 
one-dimensional optimal switching problems, where the switching times are 
restricted to lie in a discrete time set. This allows us to prove a stability 
result for obliquely RBSDEs which will be used several times in the paper. 

2.1. Definition. A discretely obliquely reflected BSDE is a reflected BSDE 
where the reflection is only allowed on a discrete time set. 

We thus consider a grid ?R. := {tq = 0, . . . = T} of the time interval 
[0,T] satisfying 

(2.1) m := max \rk - rk-i\ < -. 

l<k<K K 

We also consider a matrix valued process C = {C^-')i<ij<m such that C^^ 
belongs to 5^ for i,j£{l,...,d} and satisfies the structure condition 



(2.2) 



( Cf = 0, for 1 < i < d and < t < T; 

inf CV > — , for 1< i, j < d with i ^ j\ 

o<t<T * L - - T-J-: 

inf C'l^ + Cf - Cf > 0, for l<i,j,l<d with i^j^j^l. 
k 0<i<T * * ' - ^Ji - -r J^J T- 

We introduce a random closed convex set family associated to C: 

Qt ■■= {y G K"' I > max(y^' - C^^), 1 < i < d}, < t < T, 
and the oblique projection operator onto Qt, denoted Vt and defined by 
Vt:y€M''^(max{y^ -Ci^}) 

\ jeX / l<i<d 

which is ^ C?) ;B(M'^)-measurable. 

Remark 2.1. (i) It follows from the structure condition (2.2) that V 
is increasing with respect to the partial ordering relation y, where y ^y' 
means y* > (y')* for all i £l. 

(ii) An easy calculation leads to 

\'Pt{yi)-Vt{y2)\ < Vd\yi-y2\ for any yi,y2 GM"^. 

We observe that the constant y/d is optimal in our setting taking, for exam- 
ple, yi := (maxjj- C^-^ , 0, . . . , 0) and y2 := (maxj j Q*"' -|- 1, 0, . . . , 0) . Thus Vt is 
-L-p-Lipschitz continuous with L-p := ^/d. 

Finally, we are also given a random variable ^ G [L'^{J-t)]'^ valued in Qt, 
representing the terminal value of the BSDE and a random function F -.O, x 
[0,T] X M'^ X M"^'" M"^ which is ^ ® B(M'^) (g) B(A^'^''')-measurable and 
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satisfies the Lipschitz property 

\F{t, y, z) - F{t, y', z')\ < L{\y - y'\ + \z- z'\) 

for all {t,y,y',z,z') G [0,r] x (W^)^ x {M'^^'^)^, P-a.s. We shall also assume 
that 

(HF) The component i of F{t,y,z) depends only on the component i of 
the vector y and on the row i of the matrix z, that is, F^{t, y, z) = F^{t, y^, z^). 

Given this set of data (3ft, C,F,,^), a discretely obliquely reflected BSDE, 
denoted P(3ft, C, F, e), is a triplet (Y^, y^, Z^) G (5^ x 5^ x H'^)^ satisfying 
Y^ = Y^ ■.= ieQT, and defined in a backward manner, for j < k — 1 and 
t e [rj,rj+i), by 

Y^ = Y,^^ + F{n,Y^,Z^)du- H^'z^dWu, 



(2.3) { 7^+^ 

This rewrites equivalently for t € [0, T] as 

>f = e + i^(n, yj^, d^ - jj"" dW^ + (K^ - Kf) , 

(2.4) {Kf:= ^Kf^r<t}, 

reS\{0} 

with AKf := - = -{Yt^ - yt?). 

Observe that K"^ G (A^)^, since C*-' is nonnegative and valued in 5^, for 
any i,j el. 

We shall also use the following integrability condition for some p>2: 

(Cp) \e+ sup r|F(s,0,0)rds</3, 

te[o,T] JO 

where /3 is a positive random variable satisfying E[/3] < Cl. Importantly, /3 
does not depend on K. 

The proof of the following a priori estimates is postponed until the Appendix. 

Proposition 2.1. Assume that (Cp) holds for some given p>2, there 
exists a unique solution (y^,y^,Z^) to (2.3) and it satisfies 

2.2. Corresponding optimal switching problem. In this subsection, we in- 
terpret the solution of the discretely obliquely RBSDE (2.4) as the value 
process of a corresponding optimal switching problem, where the possible 
switching times are restricted to belong to the grid JR. Our approach relies on 
similar arguments as the one followed by Hu and Tang [15] in a framework 
with continuous reflections. 
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A switching strategy a is a nondecreasing sequence of stopping times 
(6'j)jgN, combined with a sequence of random variables (aj)jgN valued in X, 
such that Oj is -measurable, for any j G N. We denote by A the set 
of such strategies. For a = {6j,aj)ji^fq E A, we introduce iV" the (random) 
number of switches before T as 

(2.5) N"" = #{k£N* ■.0k<T}. 

To any switching strategy a = {9j,aj)j^^ G A, we associate the current state 
process (at)ie[o,T] the compound cost process (^f )tG[o,r] defined, respec- 
tively, by 

TV" TV" 

at := aol{o<t<0o} + ^aj-il{g^_^<t^g^y and A^ — X] ^? ^{«j<*<t} 
i=i i=i 

for < t < T. For {t,i) £ [0, T] x X, the set At,i of admissible strategies 
starting from i at time t is defined by 

At,i = {a = {Bj,aj)j eA\eo = t,ao = i,E[\A^\'^] < oo}. 

Similarly, we introduce A^^, the restriction to K-admissible strategies 

Af^i := {a = (0,, a,),eN G At,i \ 9j G K, Vj < iV"} 

and denote A^ := Ui<d-^o!r 

For {t,i) G [0,r] X X and a G ^f^, we consider as in [15] the associated 
one-dimensional switched BSDE defined by 

(2.6) ^" 
-A^ + Al, t<u<T. 

Theorem 3.1 in [15] interprets each component of the solution to the con- 
tinuously reflected BSDE (1.1) as the Snell envelope associated to switched 
processes of the form (2.6), where the switching strategies a are not re- 
stricted to lie in the reflection grid The next theorem is a new version 
of this Snell envelope representation adapted to the context of discretely 
obliquely reflected BSDE (2.4). 

Theorem 2.1. Assume that (C2) is in force. For any « gX and t G 
[0,r] , the following hold: 



(i) The process dominates any "^-switched BSDE, that 



IS, 



(2.7) C/f < {Y^)\ P-a.s. for any a G A\ 



■i,t- 
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(ii) Define the strategy a* = {6j,a*)j>o recursively by (^qjOq) := (t,i) 
and, for j > 1 , 

e* := inf {. G [e*^„T] n K I (yf < max {(yf )^ - c?-'}|, 

a* := minh / a*„i | (1^?)^ - C;;-' = max - C;;-'}}. 

T/ien, we have a* G ^f^j and 

(2.8) (yJR)j = C/f, P-a.s. 

(iii) The following "Snell envelope" representation holds: 

(2.9) (y^)^ = esssupf/f, P-a.s. 

Proof. Observe first that assertion (iii) is a direct consequence of (i) 
and (ii). Let us fix t G [0,T] and i £l. 

Step 1. We first prove (i). 

Set a = {Ok,ak)k>o G -^fi and the process {Y°',Z°-) defined, for s G [t,T], 

by 

V ■■= J2^y!'r"Me.<s<e,^^} + e^Ms=T} and 

fc>0 

(2.10) 

k>0 

Observe that these processes jump between the components of the discretely 
reflected BSDE (3.5) according to the strategy a, and, between two jumps, 
we have 



s 

Ok 



(2.11) 



ye\^,+ [ '^'F"=(s,y;,Z,»)ds- / Z^dWs + {Kl^^_r 
Introducing 
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for s G [t,T], and summing up (2.11) over k, we get, for t <u<T, 

y« = + f F''^ (s, Y,-, Z^) ds - r Z« dWs -A^T + K + K^- K- 

Ju J u 

Using the relation Y^ = Ve^, {Y^J for all G {0, . . . , iV}, we check that 
is increasing. Since [7° solves (2.6), we deduce by a comparison argument 
(see [21], Theorem 1.3) that C/f < Y^. Since a is arbitrary in ^f^j, we de- 
duce (2.7). 

Step 2. We now prove (ii). 
^Consider the strategy a* given above as well as the associated process 
{Y"-* ,Z"-*) defined as in (2.10). By definition of a*, we have 

which gives 

(2.12) / d{Kfrl=G and (y| j-^ - (y| )«;: + c,1";+^ = 
for ah A; G {0, ... , N"' - 1}. We deduce from (2.2) that 

Yf = + j F< (s, Yf , zf )ds- j zf dWs -Af+Af, t<u<T. 

Hence, {Y"' , Z"* ) and ([/"* , ) are solutions of the same BSDE and {Yt^Y = 
Uf. To complete the proof, we only need to check that a* G A^, that is, 
E|j4y p < oo. By definition of a* on [t,T] and the structure condition on 

the cost (2.2), we have \Af\ < maxkjti\C'f\ which gives E[|Af < Cl- 
Combining 

Af = Yf - Yf + j F^^s, Yf , zf )ds- J Zf dWs + Af 

with the Lipschitz property of F and the fact that (Y^,Z^) G (5^ x Ti'^)^ 
(recall Proposition 2.1), we get the square integrability of A^ and the proof 
is complete. □ 

Remark 2.2. Although the optimal strategy a* depends on the initial 
parameters t and i, we omit the script {t,i) for ease of notation. 



Combining the previous representation with the a priori estimates of 
Proposition 2.1 and the structure condition (2.2), we deduce the following 
estimates, whose proof is postponed until the Appendix. 
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Proposition 2.2. Assume that (Cp) holds for some given p>2, then 

rT \p/2 



E 



sup \U^ \P + 

.s6[t,T] 



<cl 



for the optimal strategy a* G A^^, {t,i) G [0, T] x X. 



2.3. Stability of obliquely reflected BSDEs. We now study the depen- 
dence on the solution with respect to the parameters of the BSDE. In the 
"abstract" setting considered, we obtain precious estimates for the analysis 
of the regularity of the solution to the discretely obliquely reflected BSDE 
as well as the convergence of the discrete-time scheme. 

We consider two discretely reflected BSDEs, with the same reflection 
grid but different parameters. For (. E {1, 2}, we consider an J^^-measurable 
random terminal condition a random L-Lipschitz continuous map (y, z) i— >• 
^F{-,y,z), satisfying (HE) and a matrix of continuous cost processes 
(^C^^)i<ij<d satisfying the structural condition (2.2). 

We suppose that the coefficients satisfy the integrability condition (C4). 
Eor e £ {1,2}, we denote by (^^y^/y^/z^) G (5^ x 5^ x n^f the solution 
of the obliquely discretely reflected BSDE V{^,^C,^F/^). 

Defining JF^ = ^Y^ - ^Y^, 6?"^ = - ^Y"^, 5Z« = - ^Z^, 6^ := 
— together with 

\6C,\^ :=max|^C7*J -2c7*%s), 



l<5F, 



max 



sup I 

y,zeR''xM'''i 



'F'\{s,y,z) 



for s € [0,T], we prove the following stability result. 

Proposition 2.3. Assume that (C4) holds. Then we have, for any t£ 

f-T 



E[\5Yf\'^]+E[\6Y;^\''] + -E 



\6Z'f\'ds 



<cAe 



\5Fs\l,ds + \5if 



+ E 



sup |(5Cr-|^ 



1/2 



Proof. The proof is divided into three steps and relies heavily on the 
reinterpretation in terms of switching problems. We first introduce a con- 
venient dominating process and then provide successively the controls on 
the 5Y^ and 5Z^ terms. 



Step 1. Introduction of an auxiliary BSDE. 

Let us define F := V ^F, ^ := V and C by C*^ := ^C'^ A "^C^K 
Observe that F satisfies (HE), C satisfies the structure condition (2.2) and 
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that (C4) holds for the data {C,F,C). We denote by {Y^,Y'^,Z^) the solu- 
tion of the discretely obliquely reflected BSDE P(5R,C, F, recalling (2.3). 

Using (HF), the definition of F and the monotonicity property of V [see 
Remark 2.1(i)], we easily obtain by a comparison argument on each interval 
[rk,rk+i), ke{0,...,K-l}, that 



(2.13) 

Recalling Theorem 2.1, we introduce the switched BSDEs associated to 
lySR^ 2y5R ^^^^ ySR denote by a = (^j, aj)j>o the optimal strategy related 
to starting from a fixed {i,t) £l x [0,r]. Therefore, we have 



2.14) (yf )^ = C/f = e'^^ + / F"»(s,C/,^1/")ds 



T 



Step 2. Stability of the Y component. 

Since a G Af^^, we deduce from Theorem 2.1(iii) that 



T 



I Ad I I Ad 



^e{i,2}, 



where A°- is the process of cumulated costs { C'^^)ij^x associated to the 
strategy a. Combining this estimate with (2.13) and (2.14), we derive 



C/f I + |/7f 



(2.15) \CYry-{'Yry\<\u^ 

Since both terms on the right-hand side of (2.15) are treated similarly, we 
focus on the first one and introduce the continuous processes T" := -|- A'^ 
and ^F"^ := ^C/" + ^A'^. Applying Ito's formula, we compute, for all t < n < T, 

rT 



(2.16) 



|-pa 



Ipd |2 I 
^ 111 



IK"- VfPds 



I pd 1-pd 1 2 
I T T I 



+ 2 



(Ft - iFt)[F'^^ (s, C/f , V/) - 'F'^^is, 'U!,'vn] ds . 
Since F = ^F V "^F and ^F is Lipschitz continuous, we also get 

< \5F s\oo + L{\Tl - ^Tl\ + \Al - ^ Al\ + \V^ - ^V^\), 0<s<T. 

Using classical arguments, we then deduce from the last inequality and (2.16) 
that 

(2.17) \r^-'Tf\''<CL(Et\\Se^\'' r\SFs\l,ds 



+ SUP EiPf-i^fp]). 

t<s<T 
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Moreover, using the inequality \x y y — y\ <\x — y\ for x,y G M and the 
convexity of the function x i— t- , we compute 



¥.t[\A1-'A1Y] 



(2.18) 



{ek<s} 



k=l 

<EJ|iV"|sup|(5C, 



2 

r-loo 



t<s<T. 



Plugging in (2.17) and recalling the definition of and F", we get 



\N^\snp\6Cr\l+ I |5F,|Lds + |5eP 
reSR it 



The exact same reasoning leads to the same estimate for \Uf — ^C/f p. There- 
fore, we deduce from (2.15) and the Cauchy-Schwarz inequality that 

E[|(2i;^)^-(iyf)^|2] 



(2.19) 



sup |5Cr|^ 



+ E 



\5F,\i,ds+m 



Using Proposition 2.2, we compute, since i is arbitrary, 
E[|2y«-iy«|2] 



(2.20) 



<Ci E 



5Fs\1ds + \5i{' 



+ E 



sup|5Cr|^ 



1/2 



S'iep 3. Stability of the Z component. 

Applying Ito's formula to the cddldg process |(5y^p and noting SK 
^K^ - '^K^, we obtain 



E 



|5y^SR|2+ r\szf\^ds+ ^ \A6K, 

" t J. ^„^rri 



5R|2 
r I 



t<r<T 



:E 



(jy^^P + 2 / (^yf (^F, ds + 2 / 6Yf d5K, 



5 s > 



where we used the fact that \6Y^\'^ - l^y^p -2(5y^(<5y*^-(5y^) = {ASK^l"^. 
Since SK is a pure jump process, we compute 



E 



rT 



< E 



a 



t<r<T,re3i 



t<r<T 



a > 0, 
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which, for a large enough and using standard arguments, leads to 

rT 



E 



t<r<T 



fnm'] + E [ r \6Fs\ids + Y \^^r\'] ) ■ 



<CL(nm^]+E 

t<r<T,re^ 

Since (2.20) holds true for any t € [0,T], we deduce 



E 



|(5Zf|2ds+ Y \^^^^ 



r I 



t<r<T 



\SFsLds 



+ E 



sup |(5Cr.|^ 



1/2 



which concludes the proof of the proposition. □ 



3. Regularity of discretely obliquely reflected BSDEs. This section is 
dedicated to the derivation of regularity properties for the solution of dis- 
cretely reflected BSDEs. These results are obtained in a Markovian diffusion 
setting. This means that the randomness of the parameter (C, F, ^) is due to 
a state process X, which is the solution of a stochastic differential equation 
(SDE). In this framework, we focus on the ^^-regularity of the compo- 
nent of the solution of the BSDEs. The main results are retrieved by means 
of kernel regularization and Malliavin differentiation arguments. Finally, we 
extend this result to the case where the diffusion X is replaced by its Euler 
scheme. 



3.1. A diffusion setting for discretely RBSDEs. Let X be the solution 
on [0, T] to the following SDE: 



(3.1) 



Xt = XQ+ [ b{Xs)ds+ [ a{Xs)dWs, 0<t<T, 
Jo Jo 



where Xq £ M"^ and {b, cj) : R"' ^ M"* x 7W™'«(M) are L-Lipschitz functions. 

Under the above assumption, the following estimates are well known (see, 
e.g., [17]): 



E 



(3.2) 



sup \Xt\P 
■t&[0,T] 



< Cl and 



sup I E 



sup 



iXf, — X,, P 



G[0,T] ^ ^u&[0,T],\u~s\<h 



1/p 
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for any p > 0. In the sequel, we shall denote by (3^ a positive random 
variable, which may change from line to line, but which depends only on 
sup^g[o,T] 1^*1 ^'^d which satisfies E[|/3''''|p] < C£ for all p> 0. Importantly, 
does not depend on K. 

Remark 3.1. Observe that, as in [1, 7] and contrary to [18], we make 
no uniform ellipticity condition on a. This allows us to treat the case of 
nonhomogenous diffusion by setting, for example, =t, t € [0,7"]. 

In this context, we are given a matrix valued maps c := (c*-') where 
^ij . _^ ^ ^j^g L-Lipschitz continuous and satisfy 



(3.3) <^ 



c"(-) = 0, forl<i<d; 
inf c*-' (x) > 0, for 1 < i,j < d with i ^ j; 

inf {c*J (x) + c^'(x) - c*'(x)} > 0, for 1 < i,j,l< d 



withi/j,j/^- 

We then introduce a family {Q{x))x&M.-m of closed convex domains: 

Q{x) := {yeW^\y'> max{y^ - d^{x))yi G l} 
I jex J 

(3.4) 

where X := {1, . . . , d}. 

We introduce the oblique projection operator 'P(x,-) onto Q(x) defined 

by 

V : {x, y) G M"' X M*^ ^ (max{y^ - (x)}) 

\ j&T J l<i<d 

Finally, we are given: 

(i) an L-Lipschitz function g:M."^ — t- such that g{x) G Q{x) for all 
x€W^, 

(ii) a generator function, that is, an L-Lipschitz map 

From now on, we shall appeal to the following assumption: 

(Hf ) the component i of f{-,y, z) depends only on the component i of the 
vector y and on the column i of the matrix z, that is, f^{-,y,z) =/*(•,?/*, z*). 

We denote by {Y^,Y^,Z^) the solution of the discretely reflected BSDE 
'D{^,c{X), f{X,-,-),g{X)) which reads on each interval [rj,rj^i), for j < k 



(3.5) 
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(3.6) 



fT 

- J Z^dWu + {K^-K!^), 



T 



5R r7'R\ 



re5R\{0} 



^{r<t} 



and 



y/ 



o<t<r, 



0<t<T. 



From (3.2), it follows that the data {c{X), f(X,-,-),g{X)) satisfies the 
integrability condition (Cp) for all p>2. We thus deduce from the proofs of 
Propositions 2.1 and 2.2, the following estimate on (Y^,Y^,Z^) and their 
associated optimal switched BSDEs, recalling Theorem 2.1. 

Proposition 3.1. There exists a unique solution (y^,y^,Z^) to (3.5) 
and it satisfies 



(3.7) 



sup I y^j 



\Zf\^ds 



p/2 



Kf\P 



</5- 



X 



yt<T. 



Moreover, for all {t,i) G [0,T] x X, the optimal strategy a* G Af^ satisfies 



(3.8) Et 



sup \U^ \P + 
5e[t,T] 



T \p/2 

Wfl^ds 



+ \A^'\P + \N'''\P 



3.2. Malliavin differentiability of {X,Y^,Y^,Z^). We shall sometimes 
use the following regularity assumption on the coefficients: 

(Hr) The coefficients b, a, g f and {c^^)ij are C^'^ in all their variables, 
with the Lipschitz constants dominated by L. 

We denote by B^'^ the set of random variables G which are differentiable 
in the Malliavin sense and such that ||G||][j'^ := ||G||^2 + /o^ 11^4^11^2 dt < 
oo, where DtG denotes the Malliavin derivative of G at time t <T. After 
possibly passing to a suitable version, an adapted process belongs to the 
subspace £a'^ ofTi^ whenever Vs G B^'^ for all s < T and ■= WVWli + 

-C<^ll^2dt < CO. For a general presentation on Malliavin calculus for 



/o 11^* 

stochastic differential equations, the reader may refer to [19]. 



Remark 3.2. Under (Hr), the solution of (3.1) is MalUavin differentiable 
and its derivative satisfies 



(3.9) 



sup|Ds-'^| 



< OO, 
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and we have 



sup\\DsXt - DsXuLp + 



(3.10) 



s<u 



sup \DtXs - DuXs 

t<s<T 



<Cl\t-u\^/^ 



for any < li < t < T. Let G G ©^'^(M'^). Since X belongs to Cl'"^ under 
(Hr) and V is Lp-Lipscliitz continuous, we deduce that V{Xt,G) £ B^'2(R'^). 
Using Lemma 5.1 in [1], we compute 

Ds{ViXt,G)y 

d 

(3.11) = '^{DsG^ - ^sCtj(^t))l{GJ-c»i(Xt)>maxf<j{G''-c«{Xi))} 



X l{GJ-c'^(Xt)>maxf>j(G«-c«(Xt))}- 

Combining (3.11), Proposition 5.3 in [10] and an induction argument, we 
obtain that (y^,y^,Z^) is Malhavin differentiable and that a version of 
{Djr^,D^Z^) is given by 

1^ 1 t 



k=l- 



(3.12) 



+ / v^r{x„{Y!'y,izfy)DuXsds 



rj+i 



+ / v,.nXs,iY!'y,{zfy)D,,{Y:^yds 



+ vj\Xs, (yf )\ (zf )*-)D„(zf )*■ d. 

for < ti < t < rj+i and j < k. Here, V^/* denotes Yle=i'^ z'- i recaU- 
ing (Hf). 

3.3. Representation of Z. For a E A^, we introduce the process A" de- 
fined by 



(3.13) 



A?,,:=exp / V,r^iXr,Y,^,Z^)dW, 



v,r'-(x,,,y,«,zf )|2 - v,r'-(x,,,y«,zf ) ]dr 



for < t < s < T. 

For later use, we remark 



(3.14) 



sup 



sup Al, <Cl 0<t<T,p>2, 

t<s<T -CP 
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and deduce from the dynamics of A that 



(3.15) 



sup 1 1 Aj,t - A? 1 1 £P + sup I AS - A^ 

a<^A^ ^ t<s<T 



CP 



u<t<T,p>2. 



Proposition 3.2. Under (Hr), there is a version of such that 



izfy = Et 



(3.16) 



+ / v,r^'(x„yf,zf)Af,A^.ds 



for {t,i) G [0,T], with a* = {Oj,Oi*)j>o G Af^ the optimal strategy given 
Theorem 2.1 and recalling (2.5). 



Proof. We fix j < k and, observing that the process a* is constant on 
the interval [6**,^^^^), we deduce from (3.12) and Ito's formula that 

A?;i)„(y,«)-.* 



KXe'{D^{Y'^r^e*^^+ I ^ V.rJ(X„n^,Zf)A«;Z)„X,d. 



for 9j <u <t < Oj^i- Combining (3.11) and the definition of a* given in 
Theorem 2.1(ii), we compute 

At;.jD^{Y''r^e;^, = Af,.^^(D„(y'^)"J+i),.^^ 

for j < K. Plugging the second equality into the first one and summing up 
over j concludes the proof. □ 

We conclude this section by providing a "weak" regularity property of 
in the general Lipschitz setting. In order to get rid of the previous assump- 
tion (Hr), we make use of kernel regularization arguments. Since this proce- 
dure is very classical, we do not detail it here precisely (see, e.g., the proofs 
of Proposition 4.2 in [7] or Proposition 3.3 in [1]). 

Proposition 3.3. There is a version of Z^ satisfying 
ft 



(3.17) 



E 



\z!}\^du 



<CL\t-s\, 



s<t<T. 



18 



J.-F. CHASSAGNEUX, R. ELIE AND I. KHARROUBI 



Proof. Combining (3.9), with (3.14), (3.16) and Doob's inequality, we 
observe that 

sup \\Zf\\^,<Cl, p>2, 
telo,T] 

holds under (Hr). Therefore (3.17) is satisfied under (Hr). As in the proof of 
Proposition 4.2 in [7], the stability results of Proposition 2.3 allow us to use 
classical Kernel regularization arguments. Since the previous estimate holds 
uniformly for the sequence of approximating regularized BSDE, the proof is 
complete. □ 



3.4. Regularity results. We consider a grid vr := {tg = 0, . . . ,t„ = T} on 
the time interval [0, T] , with modulus |7r| := maxo<j<n-i \ ti+i — ti\, such that 

KCTT. 

We want to control the following quantities, representing the 7^ ^-regularity 
of (y,Z): 



(3.18) 



E 



^7r{t)l 



and E 



\Z, 



5R 



7^ Ph/ 



where nit) := sup{ij € iT;ti< t} is defined on [0,T] as the projection to the 
closest previous grid point of vr and 



(3.19) 



1 



-E 



ti+i 



Z?ds 



i e {0, . . . ,n - 1}. 



Remark 3.3. Observe that {Zf)s<T '■= (■^^(s))s<T is interpreted as the 

best "W^-approximation of the process Z^ by adapted processes which are 
constant on each interval [tj,tj+i), for all i <n. 



Proposition 3.4. The following holds: 



-E 



T 



5R 







< sup E[|yf -y4|2]<CL|7r| 

ie[o,T] 



Proof. Observe first that 
ft 



n\Y^-Y;%n<E 



» |2l 

{t)\ 



f{XsXt,Zf)ds+ ZfdW, 

n{t) J-K{t) 



0<t<T. 



The proof is concluded combining this estimate with (3.2), Propositions 3.1 
and 3.3. □ 



We now turn to the study of the regularity of the process Z^. 



Theorem 3.1. The process satisfies 
rT 



(3.20) 



E 



\Z^ 



Zfl^ds 







<CL(|7r|V2 + ^|,r|). 
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Proof. A regularization argument as in proof of Proposition 3.3 allows 
us to work under (Hr). From Remark 3.3, it is clear that 



(3.21) 



E 



T 



\z. 



< E 


[L 







For s <T and a = {ak,Ok)k>o G Af^, £ G X, we define (l^/t),<t<T by 



N°- 

^V,c^^^,,^^iXe,)Ale,iDsX)e, 



k=l 

We now fix £ G X and denote, for u <T, by a" G the optimal strategy 

associated to the representation of (1^)^, recalling (ii) in Theorem 2.1. 
Observe that, by definition, we have 

(3.22) Ar'*' = Ar«" and a* = a", rj <t < u < rj+ij < k. 

Fix i <n, and deduce from Proposition 3.2 and (3.22) that 



(3.23) 



E[\{z!'Y-iz!^Y\'] = E[\v,i-v^:,f] 



<2(E[|l^y-l^-|2]+E[|V^--^-|2]) 

for t£ [ti,ti+i). Combining (Hr), (3.9), (3.10), (3.14), (3.15) and Cauchy- 
Schwarz inequality with the definition of V"", we deduce 

(3.24) E[\V,f - |2] < ClH'/^ ti<t< U+,,i< n. 
Since V^f^ is a martingale on we obtain 



(3.25) 



where the last inequality follows from (3.24). Combining (3.23), (3.24), (3.25) 
and summing up over i, we obtain 

(■T 



E 



\{zff-{z%) 



K-l 



< clHI' + ki \Ew^r? - - i<^/)j ■ 

Combined with (3.9) and (3.14), this concludes the proof since I is arbitrary. 
□ 
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3.5. Extension. We shall approximate the process X by its Euler 
scheme X'^ , with dynamics 

(3.26) X- = Xo + l^ + 1^ a{X:^^^)dWs, 0<t<T. 
Classically, we have the following upper-bound, uniformly in vr: 

(3.27) E 



1 i/p 

sup \X[\P <Cl, p>2. 

-0<t<T -I 



The control of the error between X and its Euler scheme X'^ is well 
understood (see, e.g., [16]) and we have 



(3.28) 



E 



sup \Xt-X[\P 

0<t<T 



1/p 



<C?KI'/', P>2. 



In this context, we denote by (y^"^ y^"^ 2'*^") the unique solution of the 
reflected BSDE P(3f?, c(X'^), /(X'^, •), 5(^'')). Our main result here is the 
counterpart of Proposition 3.4 and Theorem 3.1 when X is replaced by X^ . 



Proposition 3.5. The following hold: 

t-T 



E 







and 



E 



Jo 



\Zl 



7eu\2 



ds 



<Cl\'k\ 



<CL(|7r|i/2 + KK|) 



Proof. We only sketch the main step of the proof since it follows ex- 
actly the same arguments as the one used to obtain Proposition 3.4 and 
Theorem 3.1. 

Step 1. We use a kernel regularization argument which allows us to work 
under (Hr). In this case, we observe that X"" belongs to and satisfies 

DsX^ = a{Xl^,^) + f VMKir))DsX:^^^ dr 

J s 

+ tj2'^.<^'{x:^r))Dsx:^,^dw^ 

for s <t. One then checks (see [1], Remark 5.2, for details) that 

sup iDsX'^l < oo, 

s<T sv 
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(3.29) sup\\DsX^ - DsX^W^, 

s<u 



sup \DtX:-DuXJ\ 

t<s<T 



<Cl\t-u\'^\ 

LP 

< n < t < r. 



It is also straightforward that (y^", y^", Z*^") is Mahiavin differentiable and 
satisfies (3.12) with X'^ instead of X. 

Step 2. In order to retrieve the results of the proposition, one then follows 
exactly the same steps and arguments as the ones used in the previous 
Sections 3.3 and 3.4. □ 

4. A discrete-time approximation for discretely reflected BSDEs. We 

present here a discrete time scheme for the approximation of the solution of 
the discretely obliquely reflected BSDE (3.5). 

Recall that vr := {to = 0, ■ ■ ■ ,tn = T} is a grid on the time interval [0, T], 
such that K C vr and |7r|n < L. In the sequel, the process X is approximated 
by its Euler scheme X'^ (see Section 3.5 for details). 

4.1. A Euler scheme for discretely obliquely reflected BSDEs. We intro- 
duce a Euler-type approximation scheme for the discretely reflected BSDEs. 
Starting from the terminal condition 

y^'" = yf := g{x^) G c(XJ), 
we compute recursively, for i <n — 1, 



(4.1) 



y,«'- = ek5^ I j-^j + (t,+i -u)fixz,Y,^^\zf^n 



This kind of backward scheme has been already considered when no re- 
flection occurs (see, e.g., [4]) and in the reflected case (see, e.g., [1, 7, 18]). 
See also [3] for a recent survey on the subject. 

Combining an induction argument with the Lispchitz-continuity oi f, g 
and the projection operator, one easily checks that the above processes are 
square integrable and that the conditional expectations are well defined at 
each step of the algorithm. 

Remark 4.1. (i) This so-called "moonwalk" algorithm is given by an 
implicit formulation, and one should use a fixed point argument to compute 
explicitly Y^''^ at each grid point. 

(ii) In the two-dimensional case, Hamadene and Jeanblanc [13] interpret 
y^ — y^ as the solution of a doubly reflected BSDE. It is worth notic- 
ing that the solution of the corresponding discrete time scheme developed 
by [7] for the approximation of doubly reflected BSDE exactly coincides with 
(y!R,^)i _ (ySR,^)2 derived here. 
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For later use, we introduce the piecewise continuous time scheme asso- 
ciated to the triplet (y3?,7r y3?,7r^ ^!R,7r^_ gy ^j^g martingale representation 
theorem, there exists Z^''" € T-L^ such that 



= IE*. Kr;:] + / dWu, i<n-l 



and by the Ito's isometry, for i <n — 1, 



(4.2) 



-E 



H+l 



ti + l 



We set Zf := Z'^^^^ for t e [0,T], define F^'^ by 
(4.3) 



y«,. ^ y^j,. ^ (^^^^ _ y,f '^ zj'-) 



z^^^dWu, k<t<ti+ui£i, 



and introduce on [0,r] by := Y^^'^l^.^^y +r{X^ ,Yr''')litem- 
This can be rewritten as 



(4.4) 



y,«'- = 5(x?) + /(x;(„),y-;:),z«'-)d..- / z^dTV. 

+ (irj'"-irf'"), 

Kf'-:= Aiff'-l|,<,} and 



reK\{0} 



AKf := y^ 



5R,7r 



5R,7r 



I ^^t""'" = + P(Xf , y-'-)i{,,^|, < t < T. 

We finally provide a useful a priori estimate for the solution of the discrete 
time scheme whenever / does not depend on z, whose proof is postponed 
until Appendix A. 2. 

Proposition 4.1. /// does not depend on z and \it\L < 1, the following 
bound holds: 



(4.5) 



E 



sup \Yf. 



5R,7ri 



■0<i<n 



<Cl p>2. 



Recall that neither depends on 3? nor on ir. 

4.2. Convergence results. The next proposition provides a control on the 
error between the discrete-time scheme (4.1) and the solution of the dis- 
cretely reflected BSDE (3.5). 
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Proposition 4.2. The following holds: 



(4.6) 



sup E[|y« - + |yf - -^p] + e 

te[o,T] 



\Zl 



ds 



<C7,.|Lp|2-(|vr|V2 + ^|,,|), 

where we recall that L-p = Vd is the Lipschitz constant of the projection 
operator V. 

Proof. As in Section 3.5, we consider (y^"^ yeu^ ^eu^ ^j^g unique so- 
lution of the reflected BSDE V{?fi,c{X^)JiX'',-),g{X^)). Using Proposi- 
tion 2.3, the Lipschitz property of /, g, c and (3.28), we obtain 



1, 



(4.7) 



sup E[|yf - y/"|2 + |yf - y/"n + -e 

t€[0,T] 1^ 



T 



\z'^-ZT?ds 



<Cl\A- 

Using the same arguments as in the proof of Proposition 3.4.1, Step l.a 



in [6], for example, we get the following inequality: 

sup ] 

<e[ii,ti+i) 



sup E [| y/" - yf^^ 1 2 + 1 y/" - yf 1 2] 



+ E 



(4.8) 



z 



5R,7r|2 



ds 



<Cl E 



+ 



rU+i 
/ ^' 



y^-y3„^l' + |z^-^f(,)|')ds 



-7r(s)l 



There are two differences with the proof of Proposition 3.4.1 in [6]. First, 
V here depends both on x and y: but this is not a problem since (y^n yen^ 
Z^'^) and (y3?,7r^ y5R,7r^ ^5R,7r-j parametrized by the same forward pro- 
cess X'" . 

Second, V is not 1-Lipschitz but only Lp-Lipschitz, with L-p > 1, in its y 
component. This explains the term \L'p\'^'^ in (4.6). Indeed, we have, for 
i < n, 

This leads, using an induction argument (see, e.g.. Step l.b in the proof of 
Proposition 3.4.1 in [6]), to 

rT 



sup E[|y/" - y^^'^i^ + |y/" - y^^-'^i^] 

ie[0,T] 



+ E 



|Zr-^f'"|^ds 



24 J.-F. CHASSAGNEUX, R. ELIE AND I. KHARROUBI 

< CilLpp- f Ivrl + Tdyr - yX)\' + \ZT - Z^iJ) 







Combining the last inequality with Proposition 3.5 and (4.7) completes 
the proof. □ 

The term iLpp*^, even when k is small can be very large. Moreover, we 
shall see in the next section that it yields to a poor convergence rate for 
continuously reflected BSDEs. This term is due to the "geometric" approach, 
used in the proof of Proposition 4.2, and the fact that V is only L-p-Lipschitz 
with L-p > 1. We obtain below a better control, using the stability results 
proved at the end of Section 2 but unfortunately under the assumption 
that / does not depend on z. The optimal choice for k in terms of |7r| is 
discussed in Section 5.3 below. 



Theorem 4.1. /// does not depend on z, the following holds: 
sup ] 

t£[0,T] 



sup E[|y,« - Y^'^'l' + \Yf - Y^'^f] < ClH 



E 



< Ci(K|7r| + Itt 



l/2^ 



for \tt\ small enough. 



Proof. We use here the stability results of Proposition 2.3 setting 

(lySR^lySR^l^SR) ^ (yK^yK^^SR) ^^^^ : (s, y, ^ Ff] and {^Y^ , 

2ySR,2^jj)^(yR,.^ySR,.^^K,.)^ with ^ F : (s , y , z) ^ f (X^^^^Y^^'^]) . Com- 
bining (4.5) and Proposition 3.1 with the Lipschitz property of /, it is clear 
that (C4) holds. Applying Proposition 2.3 and (3.28), we derive, for t G [0, T], 

(4.9) 

< C, (ivrl + E|yf - Y^y^ds + I ^ - Kl/ds 

Applying the discrete version of Gronwall's lemma to estimate (4.9) rewrit- 
ten at time t = tj G tt, we deduce 

E|y« _ y«,^|2 < c, L\ + f\\Y!^ _ yJJ^^|2 d. 

(4.10) ^ * 

0<t<tj <T,tj G vr. 

Plugging this estimate into (4.9), we compute 



T 
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<Cl[ KI + 



E\Y' 



< t < T, 



which combined with Proposition 3.4 leads to the first claim of the theorem. 
Observe from the representations (3.19) and (4.2) that 



E 



+ E 









Plugging (3.20), estimate (4.9) written at time t = and the first claim of 
this theorem into this expression concludes the proof. □ 

5. Extension to the continuously reflected case. In this section, we ex- 
tend the convergence results of the scheme (4.1) to the case of continuously 
reflected BSDEs. To this end, we show that the error between discretely and 
continuously obliquely reflected BSDEs is controlled in a convenient way. 

5.1. Continuously obliquely reflected BSDEs. In the sequel, we shall use 
the following assumption on /: 

(i) (Hz) The function / is bounded in its last variable: snp ^^j^d.q\f {0,0, 
z)\ < Cl and the following assumption on the cost c. 

(ii) (He) For i,j G X, the function c*-' is equal to ^c*-^ — ^c*-' , with is 
with bounded first and second derivatives and ^c*-' is a convex function with 
bounded first derivative. 

This last assumption is needed to retrieve some regularity on the reflecting 
process K (see Lemma 5.1 below). 

We denote by {¥, Z, K) £ (5^ xTi'^ x A^)^ the solution of the continuously 
obliquely reflected BSDE C{[0,T],c{X), f(X,-),g{XT)) deflned by 



(5.1) 



Y;=g\XT)+ / f{Xs,Y:,Zl)ds- 



y/>max{y/-c^^(Xi)}, 

T 

y/ -max{y/ -c^^(Xi)} dKi = 0, i£l. 

L jeX J 



ZldWs + K'T-Kl, 
0<t<T, 



Under the assumption on f, g and c, the existence and uniqueness of such 
a solution is given in [14, 15]. 

The solution of (5.1) has also a representation property in term of switched 
BSDEs, recalling (2.6). Here, of course, the switching times of the strategy 
are not restricted to take their values in K. We refer to [8] for more details. 
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Theorem 5.1. There exists, for any fixed initial condition {t,i) € 
[0,T] X I, an optimal switching strategy a := {Ok,OLk)k>Q £ ■^t,i, such that 



(5.2) 



y/ = C/f = esssupC/^ 



"-a.s. 



We deduce from (5.2), Theorem 2.1(iii), the monotonicity property of V 
and (5.1) 

(5.3) y ^ ^ for any grid 

Moreover, most of the estimates presented in Section 2 for discretely re- 
flected BSDEs hold true for continuously reflected BSDEs. For reader's con- 
venience, we collect them in the following proposition. The proof itself is 
postponed to Appendix A.S. 

Proposition 5.1. The following a priori estimates hold. For any p>2. 



(5.4) 



T 



\Zs\^dsj 



\ P/2- 



+ Et[\KT-Ktn<Et[f3''], 



0<t<T, 



and, for all {t,i) € [0,T] x I, the optimal strategy a G At,i satisfies 
(5.5) eJ sup \U^\P] +Ei[|iV'i|P] <Ei[/3-'^]. 

5.2. Error between discretely and continuously reflected BSDEs. We first 
provide a control of the error on the grid points of K between the solutions 
of the obliquely discretely and continuously reflected BSDEs (3.6) and (5.1). 



Theorem 5.2. Under (Hz), the following holds: 



(5.6) 



E 



sup{|y,-y,»*|^ + |y,-y,«|^} <ci\r\ 

■rG3? J 



1-e 



e>0. 



Moreover, if the cost functions are constant, the last inequality holds true 
with e = 0. 

Proof. The proof of this result relies mainly on the interpretation 
in terms of switched BSDEs provided in Section 2.2. For a fixed {t,i) G 
[0,T] X I, we associate to the optimal strategy d = {Ok,ak)k S A.t,i not re- 
stricted to lie in the grid the corresponding "discretized" strategy a := 
{dk,ak)k>o e -^fi defined by 



(5.7) 



Ok := iiii{r > 9k;r £R} and ak-=ak, k>0. 



Step 1. We first derive two key controls on the distance between A"" 
and A"-. 
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We fix p > 2 and, since 9^ < 0/^, k > 1, we compute 

f T \ p/2 



^ \At-A'^fdsj 

E 



(5.8) 



<s 



k=l 



\ P/2 



ds 



fc=l 

2 X p/2 

ds I . 



ds 



T 



k=l 



Using the convexity inequality (^^=1 \xk\)^ < ^ J2k=i l^fcl^, we obtain 



T 



(5.9) 



TV" 



k=l 



2 X p/2 



ds 



<Cl{l+ sup |Xt|P)|Af<^|P|K|P/2. 



te[o,T] 

Using once again the same convexity inequahty with p = 2, the Lipschitz 
property of the maps (c*-^)jjgj and the definition of 6k and Ok, we get 



AT" 



fc=l 



TV" 



ds<c2iiv^rij;ix,,-x, 



k=l 



< Cl\N^\Px^^\'P, 



where x'^''^ := J2k=l^^Prelrk-urk] - XrjP. 



Plugging this estimate and (5.9) in (5.8), we deduce 

f-T \ p/2 

l^^-^sl'ds 



(5.10) 



<Cl\NT[[l+ sup \Xs\P 

sG[0,T] 



+ x 



»l,p 



Observe also that, for r G we have Ig = lefc<r which gives 



<r 



(5.11) 



^=1 
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Step 2. We now prove the main result of the theorem. 

We introduce the processes F :=[/" — ^4" and T := U"- — A'^ , so that 

(5.12) lU"" -U^\<\T -tl + lA" - A^\. 

Applying Ito's formula to the continuous process |r — Fp on [t,T], using 
Gronwall's lemma and the Lipschitz property of /, we obtain 

|2 



(5.13) 



n-n]{Xs,u^,v,n\'ds+ / i^^-^^i^ds 



Elevating this expression to the power |, we deduce 



(5.14) 



\[n-n]{Xs,utv,^)\'dsj 



\ P/2 



+ 



T \ p/2-| 

lA^-A^pdsJ 



Combining the definition of 9 with the Lipschitz property of / and (Hz), we 
compute 

''^\[n-nKXs,utVsn\'ds 



i: 



k=l 



ds 



<CL\N^f sup {l + \Xs\^ + \U^\^ 
se[o,T] 

Plugging the last inequality and (5.10) in (5.14), we deduce 



tt-rt\p<clEt 



\N^\P( sup {l + \Xs\P + \U^\P)\W^^ + X 
^se[o,T] 



5R|,p 



Restricting to the case where t G K, we deduce from (5.11) and (5.12) that 

|y/-(y^^)^|2 



2/p 



\'!Si\+Et[\NY\x^^h'^^^ 



\NT sup (l + |X,|P + |n|P 

sG[0,T] 

Using Cauchy-Schwarz inequality and Proposition 5.1 with the last in 
equality, we obtain 

\y; - (yf )f < C];(/3^|3f?| +/3^Et[|xl^''P|T^'')- 
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Again using Cauchy-Schwarz inequality and defining Mt := Et[|x'^''^P]! 
we get 

1/2n 



(5.15) E sup|y/- (yff|2 <Cl(\^\+K 



4eR 



sup |Mt|2/P 

■tG[0,T] 



Combining Burkholder-Davis-Gundy and convexity inequalities with (3.2), 
we compute 

e[sup iMtl^H <C£(|Mo|2/*' + E[|Mt|Y^^) <C'lIE[|x'^''^|Y^^ 
<C7£|k|^/p|K|2. 

Plugging this expression in (5.15), we deduce (5.6) from the condition < 
L and the arbitrariness of i. 

Step 3. We finally consider the particular case where the cost functions 
are constant. Following the same arguments as in Step 1, we observe that 
(5.10) turns into 

rT \p/2 

/ l^^-^^pds <CI\NYU+ sup \Xsn\w^^, 
Jt J ^ se[o,T] ' 

and that A^ — A^ = ^ for r G The same reasoning as in Step 2 then leads 
to 



\Yi -Yi\^ <C'jj^t 



\N^\^ sup (i + |x,|p + |y,nl^^^|K|. 



Using Cauchy-Schwarz and Proposition 5.1 concludes the proof. □ 

We now present the main result of this section, which allows us to con- 
trol the error between the solutions of the continuously and the discretely 
obliquely reflected BSDE at any time between and T . 

Theorem 5.3. Under (Hz)-(Hc), the following holds: 



T 



\Z,-Zf\^ds 



sup E[\Yt - Yt^f + \Yt - YtY] + E 

t€[0,T] 

<Cm\^/^-', e>0. 

//, furthermore, the cost functions are constant, the previous estimate holds 
true for e = 0. 



In order to prove this theorem, we first state the following lemma dis- 
cussing the regularity of K. 
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Lemma 5.1. Under (Hz)-(Hc), there exists some positive process ij sat- 
isfying \\rj\\-^2 < Cl and such that, for all z G X, dKl < rjsds in the sense of 
random measure. 

Proof. We follow here the main idea of the proof of Proposition 4.2 
in [10] and divide the proof in three steps. 

Step 1. Fix i,j G I. We first observe using Ito-Tanaka formula, that, 
under (He), 

c'^ {Xt) = (Xo ) + f bi^ ds + [\'J dWs - f dAi^ , < t < r, 

Jo Jo Jo 

where A*-' is an increasing process and 

(5.16) W^\\y,2 + \\U^^\^2<CL. 

We then introduce F*-' := — + d^{X) > 0. Using once again Ito- 
Tanaka formula, we compute 

[Fi^]+ = [F-i/]+ + [\-fHXs,Y:,Zl) + f\X„Yi,Zi) + b^)l^^.,^^^ds 



+ l\ul^+Zl-Zi)l^^.,^^^dWs 

+ 1^ l^^^^^i-dKl + dKi - dAl^) + dL^ 

for < i < T, where is the local time at of the continuous semi- 
martingale F*-^ . Since F*-' > and A*-' , are increasing processes, we com- 
pute 



l^ry ^0} ^ i-f{X,,Y:,Zi) + f^{X,,Yi,Zi) + bin\r^,^o} 



(5.17) 

<CLfl + |X,|+sup|y/|+ sup \bf\]ds 



+ ^{r'^=o} '^^s 

for < s < T, where we used (Hz) in order to obtain the last inequality. 

Step 2. We now prove that 
(5.18) l^p.,^ojdK^=0 

in the sense of random measure. We first observe that l|pij_Q| dKi = jl-' dKi 
with 7^ := l{r^^=o}l{y/-P^(x.,y,)=o}- ^^^eed, if l^yj _^,^^^ y^^^^^dKi were 
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a positive random measure on [0,T], this would contradict the minimahty 
condition (5.1) for K. 

Suppose the existence of a stopping time r smaller than T, such that 

(5.19) r;J =0 and yI -V\Xr,Y^) = 
By definition of the projection V, we have 

(5.20) yI - {Xr ,Yr)= - Y^^ + (9^^ [Xr), 

where kr takes value in X. Moreover, Y^ — Y^^ + c*'^^(Xt-) > 0, which leads, 
combined with (5.19) and (5.20), to c*-'(X^) + c>^-(Xr) - ^^^{Xr) < and 
then contradicts (3.3). 

Thus, = for any stopping time r smaller than T and we deduce 
that 7*-' is undistinguishable from 0, which proves (5.18). 

Step 3. To conclude, using once again the minimality condition for K 
in (5.1), observe that dKl = l|pij^Q| dKl < rjs ds, with -q := Cl{1 + \X\ + 

sup^gx|y^| + sup^fcgjl^ I) which satisfies ||ry||-^2 < Cl, recalling (3.2), (5.4) 
and (5.16). □ 

Proof of Theorem 5.3. Fix t G [0,r] and introduce 5Y ■.= Y - Y^ , 
SY ■.= Y- y^, 5Z:=Z- and 5f := /(^, Y, Z) - f{X, Y^, Z^). Applying 
Ito's formula to the cddldg process |5yp, we get 



(5.21) 



\5zyds 



\SYt\' 



6Ys- d6Ys 

(*''^] t<s<T 



Recalling that 5Ys- = ^Yg, J^^^ j.-^ 6Ys dK^ > and the Lipschitz property 
of /, standard arguments lead to 



E 



(5.22) 



T 



\6ZJ^ds 



SY,dK, 



< 



5YdK, 



Using the expression of 6Y and Lemma 5.1, we obtain 



5Ys<8Yr^^, + r^\6f^ + r]u)du- 6ZudWu, 

J S J S 



rj<s< rj+i,j < K. 
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Combining (Hz), (3.2), (3.7), (5.4) and the fact that ||?/||-^2 < Cl-, we deduce 

j<K 





<E 








'-j<K ^ 



+ E 



{5fu + nu)dudKs 



< ClM +e\ Kt sup \6Yr 



Plugging this expression in (5.22) and using Cauchy-Schwarz inequahty 
together with (5.6) and Proposition 2.1 concludes the proof. □ 

5.3. Convergence of the discrete-time scheme. Combining the previous 
results with the control of the error between the discrete-time scheme and 
the discretely obliquely reflected BSDE derived in Section 4, we obtain the 
convergence of the discrete time scheme to the solution of the continuously 
obliquely reflected BSDE. In the next theorem, we detail the corresponding 
approximation error for different optimal choices of reflection time step |K| 
with respect to the discrete time step |7r|. 



Theorem 5.4. The following hold: 
(i) // (Hf)-(Hc) holds, taking \^\ ~ 



'°f ^7 for e > 0, we have 

-elog|7r| ' 



sup E[\Yt-Y^ 

iG[0,T] 



T 



\z, - z 



ds 



< 



-log(|7r|)]V2-e- 

(ii) // / does not depend on z and |7r|L < 1, taking similar grids 3? = tt, 



we have 



supE[|y,, - '^p + \Yt^ - y-'^^n < ci\^ 

i<n 



1-e 



e>0. 



Moreover, under (He), 

sup E[\Yt - y^^-'p + \Yt - y^^'^p] < C7£|7r|V2-e, ^ > q. 
te[o,T] 

(iii) Under (He), if f does not depend on z and \7r\L < 1, taking |K| 
|7r|^/^, we get 

rT 



E 



<C7£|^|i/3-^, 



e>0. 



(iv) Furthermore, for constant cost functions, (ii) and (iii) hold true with 
e = 0. 
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Proof. For e > 0, setting such that |K| 
bining Proposition 4.2 and Theorem 5.3, that 



log L-p 
-£log|7r| ' 



we obtain, com- 



sup K[\Yt - y^^'^p + \Yt - y^'^-'^n + E 



^5R,7r|2i 



T 



ds 



-1 



1/2-e 



V TT 



1/2-e 



.logdvrl). 

Therefore, (i) is proved. Furthermore, (ii), (iii) and (iv) are direct conse- 
quences of Theorems 4.1 and 5.2 or 5.3. □ 



APPENDIX 

A.l. A priori estimates for discretely RBSDEs. We collect here the 
proofs for a priori estimates given in Propositions 2.1 and 2.2. 

Proof of Proposition 2.1. Observing that on each interval [rj , rj+i) , 
(y^, y^, Z^) solves a standard BSDE, existence and uniqueness follow from 
a concatenation procedure and [20] . The rest of the proof divides in two steps 
controlling separately and (Z^,K^). 

Step 1. Control ofY^. 

As in the proof of Theorem 2.4 in [14], we consider two nonreflected 
BSDEs bounding Y^. 

Define the M'^-valued random variable ^ and the random map F by {(^y : = 
Eti and (Py := ^ti K^)i for l<j<d. 

We then denote by {Y,Z) G (5^ x 'H'^)-^ the solution to the following 
nonreflected BSDE: 

(A.l) Yt = i + P{sX,Zs)ds- j'^ ZsdWs, 0<t<T. 

Since all the components of Y are similar, Y £C. 

We also introduce {^Y, ^Z) the solution to the BSDE 

°Yt = C+ [ F(s,0y„0z,)ds- / o^.dt^,, o<t<r. 

Jt Jt 
Using a comparison argument on each interval [rj^rjj^i) and the monotony 
property of we straightforwardly deduce ^Y ^ Y^ ^ Y . 

Since (^Y,Y) are solutions to standard nonreflected BSDEs, usual argu- 
ments lead to 

(A.2) sup |yf|*'< sup |°y,|*'+ sup |y.|^=:/3, 

0<s<T 0<s<T 0<s<T 
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where the positive random variable /3 satisfies classicahy E[/3] < C^, under 
condition (Cp) for a given p>2. 

Step 2. Control of (Z^,A'^). 

We fix t < r and applying Ito's formula to the cddlag process |y^p on 
[0, t] to derive 



\^t I 



yn^|2 + 2/ y^dyf + 



5R Av^ 



f|2_|y5|2_2y»Ayf) 



Since the last term on the right-hand side is nonnegative, we deduce that 



zf |2 ds < |yf I" + 2 / Y^iFis, yf , ) ds 



T 



+ 2 / y,*LdJC + 2 / (zfyf)dVFs- 

Using standard arguments, together with (A. 2) and (Cp) for a fixed p>2, 
we compute 

(A.3) l^f P ds < Cl + ^yP{K^ - Kf) + j\zfY!^) dW,^ . 

Moreover, we get from (2.4) and (Cp) that 

rT 



(A.4) \K^-Kf^\'^<CL 



\Zf\'ds + 



Combining (A.3) and (A.4) we obtain 
(A.5) '^^ 



'■^|Zf|2ds<^ 



Zf\Hs + e 



ZfdW, 



ZfdWs 



+ Cl / (zfyf)dw, 



for any e > 0. Elevating the previous estimate to the power p/2, it follows 
from Burkholder-Davis-Gundy inequality that 

p/2n 



Et 



T 



p/2- 
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T \ p/4-| 



sup |yf |p 



+ e 



Using (A. 2) and (Cp), we deduce, for e small enough, 



(A.6) 



T 



izfpds 



Taking (A. 4) up to the power |, and combining Burkholder-Davis-Gundy 
inequality with (A.6) yields Et[|K^ - Kf\P] < C£Et[/3], which concludes the 
proof of the proposition, recalling (Cp). □ 

Proof of Proposition 2.2. Fix {t,i) G [0,r] xX andp > 2. According 
to the identification of {U""* ,¥"•*) with {Y"-* ,Z"-* ), obtained in the proof of 
Theorem 2.1, we deduce from Proposition 2.1 the expected controls on [/" 
and V"" . Writing the equation satisfied by {U"" , ) and using standard 
arguments for BSDEs, we observe that 

r-T \ p/2-| 



M\^T r] <c£ E 



sup ic/f r + f / \yf?^s] 



By definition of a* and (2.2), we have \A1'\ < max^^j ICj'^^l, which plugged 
in the previous inequality leads to Et[|A|, \p] < CfE,t[f3], recalling (Cp). 

We finally complete the proof, noticing from (2.2) that E([|Af°'*|P] < C£ x 
Et[\A^'\P]. □ 

A. 2. A priori estimates for the Euler scheme. This paragraph provides 
the proof of Proposition 4.1, concerning a priori estimates for the Euler 
scheme associated to RBSDEs. 

Proof of Proposition 4.1. The proof follows exactly the same argu- 
ments as in Step 1 of the proof of Proposition 2.1 above. The only difficulty 
is the use of a comparison argument for Euler scheme that we provide below 
in Lemma A.l. □ 

We detail here a comparison theorem for discrete-time schemes of BSDEs 
in the case where the driver does not depend on the variable z. 

For k = l,2, let be a square integrable random variable and -0/0 -I^"* ^ 
R"^ — )• M an L-Lipschitz generator function. We suppose that > ^2 and 
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■01 — on W"' X W^. For a time grid vr, we denote by V"'^ the discrete- 
time scheme starting from the terminal condition Y^'^ := S,k and computing 
recursively, for i = n— 1,...,0, 

(A.7) y,^''= = EK:;^j + 

Lemma A.l. For any vr such that \i^\L < 1, we have Yf^'^ > ^1'^' i<n. 

Proof. Since the result holds true on the grid point tn = T and follows 
from a backward induction on vr, we just prove Ytl'^ > ^n-i' Using (A.7), 
we compute 

Ytli - = E,„_, [6 - 6 I J + (T - t„_i)An_i(r,:L\ - J) 

(A.8) 

+ A„_i, 

where A„_i := .^y,:;!^^ - M^l^X^'j > and 

A„_i := <; yr.i _ yr-s ' " ^tn-i ^tn-i ^ ^' 

I ^n— 1 ^n~l 

1 0, else. 

(A.9) 

Since tpi is L-Lipschitz, the condition \7r\L < 1, implies (T — t„_i)A„_i < 1. 
Plugging this estimate, A„_i > and > ^2 and ipi in (A.8), the proof is 
complete. □ 

A. 3. A priori estimates for continuously RBSDEs. This last paragraph 
is dedicated to the proof of Proposition 5.1. 

Proof of Proposition 5.1. The proof of (5.4) is a direct adaptation of 
the proof of Proposition 2.1. The only difference is in Step 1: we approximate 
(Y, Z, K) by a sequence of penalized BSDEs (see the proof of Theorem 2.4 
in [14] or Step 3 in the proof of Theorem 5.3) which are bounded by ^Y 
and Y. Estimate (5.5) follows from the exact same arguments as the one 
used in the proof of Proposition 2.2. □ 
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